Fluctuations of Chromodynamic Fields 
in Quark-Gluon Plasma^ 

Stanislaw Mrowczynski 

Institute of Physics, Jan Kochanowski University, ul. Swigtokrzyska 15, 25-406 Kielce, Poland 
and Saltan Institute for Nuclear Studies, ul. Hoza 69, 00-681 Warsaw, Poland 



Abstract 

Chromodynamic fluctuations in the coUisionless quark-gluon plasma are found as a solution 
of the initial value linearized problem. The stable and unstable plasmas are discussed. 
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In the quark-gluon plasma (QGP), which is on average locally colorless, chromody- 
namic fields, color charges and currents experience random fluctuations which appear to 
influence dynamics of the whole system. In the equilibrium plasmas there are character- 
istic stationary spectra of fluctuations which can be found by means of the fluctuation- 
dissipation relations. Fluctuation spectra in nonequilibrium systems evolve in time, and 
their characteristics usually depend on an initial state of the system. Our aim is to discuss 
chromodynamic fluctuations in equilibrium and nonequilibrium QGP. We are particu- 
larly interested in QGP produced at the early stage of relativistic heavy-ion collisions. 
Such a plasma is unstable with respect to chromomagnetic modes due to anisotropic mo- 
mentum distribution of quarks and gluons (partons) , see the review [1]. The instability 
growth is associated with generation of chromomagnetic fields which in turn influence 
various plasma properties. In particular, transport coefficients of such a plasma, which 
are controlled by the fluctuation spectrum of chromomagnetic fields, are then strongly 
modified |2]. 

Fluctuations can be theoretically studied by means of several methods reviewed in 
the classical monographs |3I4| . We choose the method which is clearly exposed in the 
handbook [5] . The method - applicable to both equilibrium and nonequilibrium plasmas 
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- provides the spectrum of fluctuations as a solution of the initial value (linearized) 
problem. The initial plasma state is assumed to be on average charge neutral, stationary 
and homogeneous. When the state is stable, the initial fluctuations are explicitly shown 
to exponentially decay and in the long time limit one finds a stationary spectrum of 
the fiuctuations. When the initial state is unstable, the memory of initial fluctuations 
is not lost, as the unstable modes, which are usually present in the initial fluctuation 
spectrum, exponentially grow. The fluctuations of chromodynamic fields in QGP have 
been discussed in detail in our study [6j. Here we briefiy summarize the results. 

The transport theory of weakly coupled quark-gluon plasma, which forms the basis 
of our analysis, is formulated in terms of particles and classical fields. The particles - 
quarks, antiquarks and gluons - should be understood as sufficiently hard quasiparticle 
excitations of QCD quantum fields while the classical fields are highly populated soft 
gluonic modes. The transport equation of quarks reads 



where Q(t, r, p) is the on-mss-shell quark distribution function which is NcXNc hermitean 
matrices belonging to the fundamental representation of the SU(iVc) group; the covariant 
derivative in the four-vector notation reads D^^ = — ig[A^{x), ■ ■ ■ ] and E(t, r) and 
B(i, r) are the chromoelectric and chromomagnetic fields. The symbol {...,...} denotes 
the anticommutator. Since the fluctuations of interest are assumed to be of the time 
scale, which is much shorter than that of inter-parton collisions, the collision terms are 
absent in Eq. ([1]). There are analogous transport equations for antiquark {Q{t, r, p)) and 
gluon (G(t,r,p)) distribution functions. 

The transport equations are supplemented by the Yang-Mills equations describing a 
self-consistent generation of the chromoelectric and chromomagnetic flelds by the color 
four-current j'^ = (p, j) 



where t°, with a = 1, ... , N"^ — 1 are the SU(iVc) group generators in the fundamental 
and adjoint representations. 

We consider small deviations from a stationary homogeneous state which is globally 
and locally colorless; there are no currents as well. The quark distribution function of 
this state is Q^m(p) — n{\)) S™\ Due to the absence of color charges and currents in 
the stationary and homogeneous state, the chromoelectric E(t, r) and chromomagnetic 
B(<, r) fields are expected to vanish while the potentials A°{t, r),A{t, r) are of pure gauge 
only. Since the plasma under considerations is assumed to be weakly coupled with the 
perturbative vacuum state, the potentials can be gauge away to vanish. 

We write down the quark distribution function as Q{t, r, p) = Q°(p) + SQ{t, r, p), and 
we assume that ^ \6Q\ and jVpQ'^l ^ |Vp(5(5| with the analogous formulas for 
antiquarks and gluons. We linearize the transport ([T]) and Yang-Mills equations in the 
deviations from the stationary homogeneous state. We assume that 6Q, E, B, A'^ and A 
are all of the same order. The linearized transport equation is 



(DO + V B)Q{t, r, p) - |{E(i, r) + v x B{t, r), VpQ{t, r, p)} = 
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The Yang-Mills equations get after the linearization the familiar form of Maxwell equa- 
tions of multi-component electrodynamics. 

The linearized transport and Maxwell equations are solved with the initial conditions 
SQ{t = 0,r,p) = (5go(r,p), E(i = 0,r) = Eo(r), B{t = 0, r) = Bo(r), by means of the 
one-sided Fourier transformation defined as 

CXD 

/(^,k) = JdtJ dVe^("*-''--)/(t,r) . 



The chromoelectric field, which solves the equations, is found as 

/ "l aJV„'(k.p)+i^Ei.(k)-i(kxB.o(k))' , (2) 

J (2ny ix) — k • V 

where /(p) ee n(p) + n(p) + 2iV,ng(p) and 6NS{r, p) = Tr [r"^ (,5go (r, p) - 5Qo{r, p)) + 
T°i5Go(r, p)] ; e*-' (a;,k) is the chromodielectric tensor of, in general, anisotropic plasma 
in the coUisionless limit; £*-'(a;, k) does not carry any color indices as it corresponds to a 
colorless state of the plasma. 

When the plasma stationary state is isotropic, the dielectric tensor can be expressed 
through its longitudinal (£i(ijj, k)) and transverse (eL(w, k)) components and the matrix 
E*-3' (w, k) = -k'^d'^ + k'k^ + Lo'^e'^ (cj, k) from the left-hand-side of Eq. Q can be inverted 
as 

(E-i)'^(c.,k) = , ) .. ^ + , .J S'^ ^) . 

cj'^eL^a;, k) k-^ w^^etIw, k) - k^ V k"' / 

Then, Eq. ^ provides an explicit expression of the chromoelectric field. Using the 
Maxwell equations, the chromomagnetic field, color current and color density can be 
all expressed through the chromoelectric field. 

The correlation functions (S* (ti, ri)£^^ (t2, r2)), (-B^(ii, ri)i?^ (^2, r2)}, where (• • •) de- 
notes averaging over statistical ensemble, are determined by the initial correlations such 
as (<57V-(ri,pi)<57Vo^(r2,P2)), {Eloin)El^{r2)), ((57V-(ri, pi)£;^(r2)) which can be aU ex- 
pressed, using the Maxwell equations, through the correlation function of the distribution 
functions. The latter one is identified with the respective correlation function of the clas- 
sical system of free quarks, antiquarks and gluons which on average is stationary and 
homogeneous. For quarks the correlation function is 

(5Q""(ti, n, pi)5QP^(i2, r2, P2))free - (5™'-<5"P (2^)3^(3) _ 

X(5(3)(r2-ri-vi(i2-ti))n(pi). 

In the case of equilibrium plasma where all collective modes are damped, we consider 
the times which arc much longer than the decay time of collective excitations. Then, the 
correlation function of the chromoelectric fields equals 
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where the fluctuation spectrum is 



{KEi) 



k2 |w2£i(u;,k)|2 V k2 / |cj2ey(a,.,k) -k2|2_ 

As seen, the fluctuation spectrum has strong peaks corresponding to the collective modes. 

As an example of a nonequilibrium system, we discuss fluctuations of longitudinal 
chromoelectric fields in the two-stream system which is unstable with rc;spcct to longi- 
tudinal modes. Nonequlibrium calculations are much more diflicult than the equilibrium 
ones. The flrst problem is to invert the matrix E'-^(u;,k). In the case of longitudinal 
electric field, which is discussed here, the matrix is replaced by the scalar function. The 
distribution function of the two-stream system is chosen to be 

/(p) = {2nfn [(5(3) (p - q) + 5^^) (p + q) 

where n is the effective parton density in a single stream. There are four roots ±uj± (k) of 
the dispersion equation £l(w, k) = 0. The solution w+(k) represents the stable modes and 
w_(k) corresponds to the well-known two-stream electrostatic instability for k2(k-u)2 < 
2y(i2(k2 — (k ■ u)2) where u = q/|q| is the stream velocity and fi"^ = gf^n/(2|q|). Then, 
uj- (k) = with < 7k G -R. 

The correlation function of longitudinal chromoelectric fields generated by the unstable 
modes is found as 



{El{tl, T^l)El{t2, r2))unstable 
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d^k e^Mri-r2) (72 + (k.„)2y 



(27r)3 



k2 



(a;^-a;2)2^2 



(7^ + (k • uf) cosh (7k(ti + t2)) + (7k - (k ■ u)2) cosh {'juih - t2)) 

As seen, the correlation function of the unstable system is invariant with respect to 
space translations - it depends on the difference (ri — r2) only. The plasma state, which 
is initially on average homogeneous, remains like this in course of the system's temporal 
evolution. The time dependence of the correlation function is very different from the 
space dependence. The electric fields exponentially grow and so does the correlation 
function both in (ti -1-^2) and (ti — ^2)- The fluctuation spectrum also evolves in time, as 
the growth rate of unstable modes is wave- vector dependent and after a sufficiently long 
time the fluctuation spectrum is dominated by the fastest growing modes. 
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